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ABSTRACT 

■ The importance of quantum effects for exotic nuclear shapes is demonstrated. Based 

t on the example of a sheet of nuclear matter of infinite lateral dimensions but finite 

thickness, it is shown that the quantization of states in momentum space, resulting from 
[ the confinement of the nucleonic motion in the conjugate geometrical space, generates 

^""^ ■ a strong resistance against such a confinement and generates restoring forces driving 

the systems toward compact geometries. In the liquid drop model, these quantum 
effects arc implicitly included in the surface energy term, via a choice of interaction 
parameters, an approximation that has been found valid for compact shapes, but has 
not yet been scrutinized for exotic shapes. 

> 

In recent years, noncompact nuclear geometries of bubbles, tori, and sheets have 
attracted considerable interestSSi'Bi'iSi in the context of nuclear multifragmentation 
studies. According to the scenarios considered,iiffii0i it has been suggested that 
nuclear systems may assume transiently exotic shapes, and then undergo a charac- 
teristic multifragment decay. One of the prominent cases of noncompact geometries 
is that of an infinite sheet. This case affords a high degree of computational sim- 
plicity in theoretical modeling attempts but shows enough features common to many 
exotic geometries to serve as a test ground for the validity of various concepts. As 
an example, it has been claimedi'S that sufficiently thin sheets of nuclear matter, 
formed dynamically during a heavy-ion collision, are subject to a new form of insta- 
bility driven by the proximity interaction of the opposing surfaces. More recently, the 
concept of this sheet instability has been applied^ to assess the stability of Coulomb 
bubbles in general and against a crispation mode, in particular. In both these latter 
cases, ifli, the analysis relies critically on the liquid-drop model (LDM). However, 
the assumptions of the LDM and, in particular, the various proposed sets of model 
parametersi00'lii were shown to be valid only for regular nuclear geometries but 
not for exotic ones. Similarly, the or Landau Vlasov method. E use cflilS 

in most of the theoretical discussions of noncompact geometries are not designed to 
handle quantum effects resulting from strong spatial constraints associated with such 
geometries. 

Below, the importance of quantum effects for the properties of nuclear matter in 
noncompact spatial distributions is demonstrated for the case of an infinite sheet. 
The scope of this study is limited to the calculation of the volume energy of symmet- 
ric nuclear matter. It is clear, however, that also the surface energy would be affected 



by quantum effects of the type considered, and that the surface energy needs to be 
calculated accordingly, before realistic model predictions can be made for the geome- 
tries of interest. Also, for the sake of simplicity, the Coulomb energy is disregarded 
in the present study. Central to the approach used in this work is the notion of bulk 
matter as opposed to surface matter. The former is characterized by a spatial uni- 
formity in the controlling parameters, and most notably, in the nucleonic momentum 
distribution, while the latter is characterized by a nucleonic momentum distribution 
that is changing rapidly with spatial coordinate perpendicular to the surface. Only 
bulk matter is the subject of the present study. 

The energy per nucleon of bulk matter can be calculated using a formalism sim- 
ilar to that suggested by Seyler and Blanchardi and employed successfully in the 
development of the droplet model.00ll! This formalism, based on the Thomas- Fermi 
approximation, was modified here to account approximately for effects of a quanti- 
zation of the nucleonic momentum component p z perpendicular to the x-y surface 
plane of the sheet. This quantization is a necessary consequence of the spatial con- 
finement of the nucleonic motion by the geometry of a sheet of a finite thickness d. 
As discussed further below, it is necessary to distinguish between the model thickness 
d, used here to construct the momentum distribution of nucleons, and the physical 
matter thickness d m describing the profile of the corresponding spatial distribution of 
nuclear matter. The former quantity d represents the width of an idealized (square) 
confining potential well with infinitely high walls and defines the finite elementary 
quantum of the perpendicular momentum component p z 

Ap, = A . (1) 

As a result of the above quantization of p z , single-nucleon states for bulk sheet mat- 
ter populate in nucleonic momentum space, discrete, infinitely thin sheets at discrete 
values of the perpendicular momentum component, p z = kAp z , where k is any nonzero 
integer. This type of population is in a clear contrast with the uniform population of 
the Fermi sphere that is usually assumed in Thomas-Fermi calculations!00lll mod- 
eling nuclear matter. An example of such discrete population of the Fermi sphere of 
a radius equal to the Fermi momentum pp is depicted in Fig. 1. 

The quantization of the perpendicular component of the momentum can be be 
described by introducing a population function /(p, z), defined as the density of 
nucleonic states in momentum space per unit nuclear volume: 



/(p, z) = ^r££i*(p, ± \jk^Pl - s{z) ) . (2) 

Here S() denotes Dirac delta function, the factor 4 represents the spin-isospin degener- 
acy, and the function S(z) describes the z-dependence of the momentum-independent 
part of the effective single-nucleon potential 
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S(z) = mod(U(z)-U(0),^), (3) 

where M denotes the average nucleon mass (taken as M=938.903 MeV/c 2 ). The 
origin of the z coordinate is set half way between the sheet surfaces. A boldface font 
is used in Eq. |2| and throughout this paper to denote a vector quantity. 




Fig. 1. Population of nucleonic momentum space for bulk matter of an idealized sheet, 
infinite in x — y dimensions. The z-axis is perpendicular to the surface of the sheet. 



For the nuclear matter density p{z), the single-particle potential U(p ; ^j, and the 
energy per nucleon ey(z), one writes then in a close analogy to Refs. i00El 



p(z) = / dpf(p,z) , (4) 

JFS(z) 

U(p,z)= fdv'f dp7(p',*M|r-r'|,|p-p'|), (5) 

J JFS(z) 



and 



ev(z) = 4r /„, dpf(p,z)[^- + lu(p,z)] , (6) 
p\z) Jfs(z) 2M 2 

respectively. In the above equations, V(r,p) denotes the nucleon- nucleon interaction, 
and / dr and Jps(z) denote integration over the (infinite) nuclear volume and over 
the Fermi sphere of a radius Pf{ z ) in momentum space, respectively. Assuming 
that the range of the nucleon-nucleon interaction is small as compared to the linear 
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dimensions of the bulk domain, the spatial integration in Eq. [5] is effectively limited 
to this domain. Then for the bulk matter of interest here, all quantities in Eqs. 0-|6| 
are z- independent, in agreement with the definition of the bulk matter. 

The nucleon-nucleon interaction was assumed to be momentum dependent and 
given by the Seyler-Blanchardcl formula 

p -r/a 2 

y(r,P) = -C T —Al-^), (7) 
(r/a) o z 

where C represents the strength of the interaction, the parameter a represents the 
range of the Yukawa force, and b denotes a critical value of the relative momentum, 
beyond which the force becomes repulsive. The quantities r and p are the distance 
between the nucleons and their relative momentum, respectively. Values of the param- 
eters of the interaction were taken to be equal to0 C = 328.61 MeV, a = 0.62567 fm, 
and b = 392.48 MeV/c. For infinite symmetric nuclear matterS, these values assure 
a volume energy per nucleon of ey=-15.677 MeV and a kinetic Fermi energy of 33.138 
MeV. 

A straightforward analytical integration in Eqs. f|, [5|, and |], using Eq yields for 
the bulk matter (the z argument is left out for the sake of brevity) 

P = j^(PzfP 2 f - ^pIf - ^Ap z p 2 zF ) - ^(Ap z ) 2 p zF , (8) 

2 

U(j,) = V + V 1 ^, (9) 

and 

1 M 

Here, 



is the average kinetic energy of a nucleon and 

b 2 

M eff = M— (12) 

is the effective mass summarizing effects of the interaction component quadratic in 
nucleonic momentum p. The strengths of the momentum-dependent and momentum- 
independent components of the effective single-nucleon potential, respectively, are 
given by 

V l = A7iCa 3 p (13) 
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and 



2M x 

V = -Vi{l--p-er). (14) 

The quantity p Z F in Eqs. |8] and [TT| denotes the maximum value of the perpendicular 
momentum p z allowed for a given Fermi momentum pp 

PzF =Pf~ mod(p F , Ap z ) . (15) 

Note that the corresponding equations for infinite nuclear matter, characterized 
by a uniform population function f(p,z) = 4/h 3 , can be readily obtained from the 
above equations || - [13] by setting p z f=Pf, while dropping all terms containing powers 
of Ap z . 

The quantity of interest in the present study is the minimum value of the energy 
per nucleon, ey, for bulk nuclear matter confined to the geometry of a sheet of finite 
thickness. This value can be obtained by varying the input value of the Fermi mo- 
mentum, pf, in a search routine minimizing ey. Note that a calculation of ey for 
given p F and d entails the use of equations [l], [l|, || [13|, [11], [12], and [10| in an ordered 



sequence. Note also that in the above equations, momenta are expressed in units 
of MeV/c, energies in units of MeV, lengths in units of fm, and masses in units of 
MeV/c 2 . Accordingly, the Planck constant is in units of MeV fm/c, /i=1239.86 MeV 
fm/c. 

Results of the calculations are summarized in Fig. 2, where values of selected 
parameters characterizing properties of the bulk sheet matter are plotted versus the 
matter thickness, d m , of the sheet. The latter quantity was evaluated for any given 
model thickness d based on the idealized matter density profile 7(2) along the direction 
perpendicular to the sheet surface. The procedure of evaluating d m is illustrated in 
Fig. 3. The density profile of interest was generated by (weighted) summing of the 
sin 2 [27i(z — d/2)p z /h] functions for the actual distribution of p z . The quantity d m was 
defined via the requirement that the density r y(d m /2) be equal to one-half of the bulk 
density where the latter was defined via the "outermost" (i.e., highest in z < d/2) 
solution of the equation 

%(z)dz = Z 1 (Z) =Z lb . (16) 





The difference between the model and matter thicknesses calculated in the above 
manner depends on d. It increases linearly from to its maximum value of approx- 
imately 2.8 fm, with d increasing from to 5.6 fm. For d > 5.6 fm, this difference 
decreases quasi-hyperbolically with thickness, reaching saturation at approximately 
2 fm for d > 20 fm. 

The main result of the present study is displayed in the top panel of Fig. 2. 
As seen in this panel, the maximum possible binding energy per nucleon decreases 
dramatically as the thickness d m of the sheet decreases. This energy becomes negative, 
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i.e., ey becomes positive and the system becomes unbound for thicknesses of about 
d m ~ 2.8 fm and below (d < 5.6 fm). This effect of an increase in energy per nucleon 
with decreasing thickness of the sheet results purely from the quantization of the 
perpendicular component p z of the momentum (see Eq. |2|). 
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Fig. 2. From top to bottom: volume energy per nucleon, compressibility coeffi- 
cient, reduced mass, nucleon kinetic Fermi energy, matter density, and strength of 
the momentum-independent single-particle potential, as functions of the sheet thick- 
ness d m . The solid curves are obtained with the Seyler-Blanchard 8 nucleon- nucleon 
interaction parameterization with values of the parameters as indicated in the legend. 
Results obtained with the parameterization of the interaction as proposed in Ref. 11 
are shown with dotted lines. 
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This feature is only approximately accounted for by the standard droplet-model for- 
mula relying on a uniform, Thomas- Fermi population of the Fermi sphere in momen- 
tum space. The latter model "adjusts" the parameters of the effective nucleon-nucleon 
interaction so as to have the average shell effects included in the surface energy term. 
The model relies here on the fact that for compact shapes (e.g., of a square box) the 
integrated (over the nuclear volume) shell effect is to a good approximation propor- 
tional to the surface area. However, the large magnitude of the volume shell effect 
raises the question whether the standard droplet model accounts accurately for these 
effects in the casesof exotic shapes and whether it is, therefore, applicable at all to 
nuclear matter in noncompact geometries. The present paper does not answer fully 
this question. At any rate, this quantum effect generates a strong resistance of nuclear 
matter against the development of noncompact geometries associated with exotic spa- 
tial confinements of the nucleons motion, in addition to that generated by a "true" 
surface tension. It gives rise to strong effective forces, akin to repulsive adiabatic 
forcesjl! driving the nuclear system away from a formation of noncompact shapes 
attained dynamically in nuclear reactions. 



1.5 
1.0 
0.5 



0.0 

1 2 3 4 5 6 

z (fm) 

Fig. 3. The relationship between the model and the matter thickness of a sheet of 
nuclear matter. 

From the top panel of Fig. 2, one notes that the value of the energy per nucleon 
for the bulk sheet matter differs significantly from the asymptotic value of eyoo=-15.7 
MeV already for sheet thicknesses d m comparable in magnitude to nuclear diameters. 
This indicates that the parameters of the Seyler-Blanchard interaction§0 are not well 
suited for absolute, more accurate calculations of the type reported in the present 
work. The same is true with respect to any parameterization that is set up so as to 
describe average trends in nuclear masses and, hence, accounts for the average shell 
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effects due to a particular shape of a sphere, but not those associated with arbitrary 
shapes. 

For comparison, in Fig. 2 results obtained with a recently proposed^'" param- 
eterization of the nucleon-nucleon interaction are shown as dotted lines. This new 
parameterization adds two terms to the Seyler-Blanchard form — an attractive inter- 
action inversely proportional to the relative momentum, and a repulsive component 
proportional to the nuclear matter density to the power of two-thirds. The integrals 
involving the former term were evaluated numerically. It is clear from the top panel 
of Fig. 2 that the quantum effect discussed above depends on the nucleon-nucleon 
interaction but that it is strong for eitherSBS parameterization. 

The remaining four panels of Fig. 2, shown mostly for the sake of completeness, 
illustrate effects of a quantization of the perpendicular momentum p z on selected 
properties of nuclear matter at the stationary density minimizing the energy per 
nucleon of bulk matter. Sharp dips and peaks in the respective functions result from 
the discontinuities in first derivatives of the underlying momentum distribution at 
values of momenta p that are integer multiples of the the elementary quantum Ap z . 
One notes, that the new parameterizationlldta results in asymptotic values of the 
parameters that are different from those obtained using the "old"H2l parameterization, 
in agreement with those reported in Ref. El In particular, the values of effective masses 
obtained with the new parameterization appear much more realistic. 

In conclusion, the present analysis demonstrates the importance of quantum ef- 
fects for noncompact nuclear shapes. These effects are shown to generate resistance 
of the nuclear systems against the development noncompact geometries and to gen- 
erate forces driving these systems toward compact shapes. The large magnitude of 
the effects discussed above allows one to question the validity of the approximation 
made in the standard liquied drop model, in which the average (volume) shell effects 
of the type considered here are implicitly included in the surface energy term. It 
allows one also to question the applicability of BUU type of computations in cases 
where noncompact geometries are involved. One notes in this latter respect that, 
similarly to the Thomas-Fermi method discussed in this study, the BUU equations 
do not consider effects of spatial confinement on the spectrum of allowed states in 
momentum space. While it is not clear, to what extent semiclassical models such as 
the modified Thomas-Fermi approach followed in the present work, can capture the 
essential features of strongly quantized systems, more accurate calculations, including 
surface and Coulomb energies, as well as effects of a finite nuclear temperature seem 
desirable. A thorough investigation of the implications of the findings made in the 
present work for several issues, appear clearly warranted. Problems of high relevance 
to current reaction studies include question as to the validity of BUU calculations, 
as well as issues associated with the magnitude of realistic droplet-model parameters 
for nuclear matter in general and not only for the specific case of spherical geometries 
of real nuclei. For example, of relevance in this context are the issues concerning 
the relative stability of the neck matter between the interacting nuclei and the decay 
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modes of this matter. This problem is especially interesting for an understanding of 
the dynamical production of intermediate-mass fragments^ in heavy-ion reactions. 

Illuminating discussions with Dr. W.J. Swiatecki are gratefully acknowledged. 
This work was supported by the U.S. Department of Energy grant No. DE-FG02- 
88ER40414. 



1. References 

[1] C.Y. Wong, Ann. Phys. (N.Y.) 77, 279 (1973). 

[2] W. Bauer, G.F. Bertsch, and H. Schultz, Phys. Rev. Lett. 69, 1888 (1992). 
[3] L.G. Moretto, Kin Tso, N. Colonna, and G.J. Wozniak, Phys. Rev. Lett. 69, 
1884 (1992). 

[4] L.G. Moretto and G.J. Wozniak, Annu. Rev. Nucl. Part. Sci. 43 379 (1993), and 
references therein. 

[5] B. Borderie, B. Remaud, M.F. Rivet, and F. Sebille, Phys. Lett. B 302, 15 

(1993) ; ibid B 307, 404E (1993). 

[6] G. Batko and J. Randrup, Nucl. Phys. A 563, 97 (1993). 

[7] H.M. Xu C.A. Gagliardi, R.E. Tribble, and C.Y. Wong, Phys. Rev. C 49, 1778 

(1994) . 

[8] L.G. Moretto, K. Tso, and G.J. Wozniak, Phys. Rev. Lett 78, 824 (1997). 

[9] R.G. Seyler and C.H. Blanchard, Phys. Rev. 124, 227 (1961); 131, 355 (1963). 
[10] W.D. Myers and W.J. Swiatecki, Ann. Phys. 55, 395 (1969). 
[11] W.D. Myers and W.J. Swiatecki, Ann. of Phys. 204, 401 (1990); ibid, 211, 292 
(1991). 

[12] W.D. Myers and W.J. Swiatecki, Nucl. Phys. A 601, 141 (1996). 

[13] G.F. Bertsch, H. Kruse, and S. Das Gupta, Phys. Rev. C 29, 673 (1984); G.F. 

Bertsch and S. Das Gupta, Phys. Rep. 160, 189 (1988). 
[14] C. Gregoire, et al, Nucl. Phys. A 465, 317 (1987). 
[15] W. Ndrenberg, Phys. Lett. B 104, 107 (1981). 
[16] J. T5ke, et al, Phys. Rev. Lett. 77, 3514 (1996). 



9 



